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Abstract 

In this paper we proved that the sequence generated by the proximal point method, associ- 
ated to a unconstrained optimization problem in the Riemannian context, has finite termination 
when the objective function has a weak sharp minima on the solution set of the problem. 
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1 Introduction 

Consider the following minimization problem 

(P) min/O) 

s.t. p e M, { ' 

where M is a complete Riemannian manifold and / : M — > M is a function. For a starting point 
p° £ M, the exact proximal point method to solve optimization problems of the form ([1]) generates 
a sequence {p k } C M as follows: 



e argmin peA/ {/(p) + \ k d 2 (p,p k )} , (2) 



where {Xk} is a sequence of positive numbers and d is the Riemannian distance (see Section [2] 
for a definition). This method was first considered in this context by Ferreira and Oliveira [I], 
when M is a Hadamard manifold (see Section [2] for a definition) and / is convex. They proved 
that, for each k € N, the function / + d 2 (.p k ) : M — » R is 1-coercive and, consequently, that 
the sequence {p k } is well-defined, with p k+l being uniquely determined. Moreover, assuming that 
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SiSV^fc = +°° an d that / has a minimizer, the authors proved that the sequence {f(x k )} 
converges to the minimum value and the sequence {x k } converges to a minimizer point. Li et 
al. [2] extended this method for finding singularity of a multivalued vector field and proved that 
the generated sequence is well-defined and converges to a singularity of a maximal monotone vector 
field, whenever it exists. 

This paper is part of a wider research program consisting of the extension of concepts and 
techniques of the Mathematical Programming of the Euclidean space M. n to Riemannian manifolds. 
It is noteworthy that in a number of recent research papers, several ideas, techniques and algorithms 
of Euclidean spaces have been extended to Riemannian manifolds and have been used for both 
theoretically and practical purposes; see [H[3l[ll[5l[6l[7l[8l[9l [TOl HH H21 [13 HH US HQ HZl [HI 
UHl [2Ql Ell 1221 (231 [2H EH] and the references therein. We observe that these extensions allow the 
solution of some nonconvex constrained problems in Euclidean space. More precisely, nonconvex 
problems in the classic sense may become convex with the introduction of an adequate Riemannian 
metric on the manifold (see, for example [18] ) . 

Following the ideas of Ferris [26J, we proved in this paper that the sequence generated by the 
proximal point method associated to the problem ([1]) has finite termination when the objective 
function is convex and the solution set of the Problem Q] is a set of weak sharp minimizers for /. 
As far as we know, the notion of sharp minimizer was introduced by Polyak [27] for the case of 
finite-dimensional Euclidean spaces; see also [28, page 205]. In this particular case it is know that 
a necessary and sufficient condition for p be sharp minimum is that € intd f(p). Rockafellar [29] 
showed that, in a space with linear structure (Hilbert space), this is a sufficient condition for finite 
termination of the proximal point method. Afterwards, Burke and Ferris [30J extended the notion 
of sharp minima to what became known as weak sharp minima, mainly to include the possibility of 
multiple solutions, and extended the previous necessary and sufficient condition for characterize the 
solution set of a minimization problem as a set of weak sharp minimizers. Li et al. [25] extended 
the notion of weak sharp minimizer to optimization problems on Riemannian manifolds as well as 
the previous result which relates finite termination of the proximal point method with weak sharp 
minima, summarized in Proposition 13.21 

The organization of our paper is as follows. In Section [21 some notations and results of Rie- 
mannian geometry as well as some fundamental properties and notations of convex analysis on 
Hadamard manifolds, are presented. In Section [31 it is presented the definition of weak sharp min- 
ima as well as some of basic related results, and proved the main resulted of the paper. Finally, in 
Section HI we made some last considerations. 

2 Notation and terminology 

In this section we introduce some fundamental properties and notations on Riemannian geometry 
and convex analysis on Hadamard manifolds which will be used later. 
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2.1 Preliminaries on Riemannian Geometry 

In this section we introduce some fundamental properties and notations on Riemannian geometry. 
These basics facts can be found in any introductory book on Riemannian geometry, such as in [31] 
and [32]. 

Let M be a n-dimentional connected manifold. We denote by T p M the n-dimentional tangent 
space of M at p, by TM = U p& mT p M tangent bundle of M and by X(M) the space of smooth vector 
fields over M. When M is endowed with a Riemannian metric ( , ), with the corresponding norm 
denoted by || ||, then M is now a Riemannian manifold. We denote by B p := {v € T p M : \\v\\ < 1} 
the closed unit ball of T p M. Recall that the metric can be used to define the lenght of piecewise 
smooth curves 7 : [a, b] — > M joining p to q, i.e., such that 7(a) = p and 7(6) = q, by 

1(7)= f b \W(t)\\dt, 

J a 

and, moreover, by minimizing this length functional over the set of all such curves, we obtain a 
Riemannian distance d{p, q) which induces the original topology on M. Given a nonempty set 
D C M, the distance function associated with D is given by 

M 3 p 1 — > d D {p) := M{d(q,p) : q £ M} € M + . 

The metric induces a map / h-> grad/ € X(M) which associates to each smooth function on M 
its gradient via the rule (grad/, X) = df(X), X 6 X{M). Let V be the Levi-Civita connection 
associated to (M, ( , )). A vector field V along 7 is said to be parallel if Vy V = 0. If 7' itself is 
parallel we say that 7 is a geodesic. Given that geodesic equation V y 7' = is a second order 
nonlinear ordinary differential equation, then geodesic 7 = 7«(-,p) is determined by its position 
p and velocity v at p. It is easy to check that j^'H is constant. We say that 7 is normalized if 
1 1 7' 1 1 = 1. The restriction of a geodesic to a closed bounded interval is called a geodesic segment. 
A geodesic segment joining p to q in M is said to be minimal if its length equals d(p, q) and this 
geodesic is called a minimizing geodesic. A Riemannian manifold is complete if geodesies are defined 
for any values of t. Hopf-Rinow's theorem ([32j Theorem 1.1, page 84]) asserts that if this is the 
case then any pair of points, say p and q, in M can be joined by a (not necessarily unique) minimal 
geodesic segment. Moreover, (M, d) is a complete metric space and bounded and closed subsets are 
compact. Take p G M. The exponential map exp p : T p M — > M is defined by exp p v = j v (l,p). 

Theorem 2.1. Let M be a complete, simply connected Riemannian manifold with nonpositive 
sectional curvature. Then M is diffeomorphic to the Euclidean space M n , n = dimM. More 
precisely, at any point p € M , the exponential mapping exp p : T p M — > M is a diffeomorphism. 

Proof. See Lemma 3.2 of [31J, p. 149 or Theorem 4.1 of [32], p. 221. □ 

A complete simply connected Riemannian manifold of nonpositive sectional curvature is called 
a Hadamard manifold. The Theorem 12.11 says that if M is Hadamard manifold, then M has the 
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same topology and differential structure of the Euclidean space R n . Furthermore, some similar 
geometrical properties of the Euclidean space R n are known, such as, given two points there exists 
an unique geodesic that joins them. In this paper, all manifolds M are assumed to be Hadamard 
finite dimensional. 

Take p € M. Let exp~ l : M T p M be the inverse of the exponential map which is also C°°. 
Note that d(q , p) = 1 1 exp~ 1 5 1 1 , the map d 2 ( . ,p) : M — > R is C°° and 

graded 2 (q,p) = -exp^p, 
see, for example, Proposition 4.8 of [32] page 108. 
2.2 Convexity on Hadamard manifold 

In this section, we introduce some fundamental properties and notations of convex analysis on 
Hadamard manifolds that will be used later. These properties can be found, for instance, in 

[331 EH ESI E]- 

A set f] C M is said to be convex if any geodesic segment with end points in Q is contained 
in £1. A function / : M — > R is said to be convex if for any geodesic segment 7 : [a, b] — > M the 
composition / o 7 : [a,b] — > R is convex. Take p € M. A vector s S T p M is said to be a subgradient 
of / at p if 

f(q)>f(p) + (s,e^v p 1 q), 

for any q € M. The set of all subgradients of / at p, denoted by df(p), is called the subdifferential 
of / at p. It is known that if / is convex then df(p) is a set non-empty, convex and compact, for 
each p £ M. In particular, given p,q £ M and u € df(p), v £ df(q), we have 

(u^xpp 1 ^) <f(q)-f(p) < (Uj-expg 1 ^). (3) 

But this tell us, from next definition, that if / is convex then df is a monotone vector fields on M. 

Definition 2.1. Let Q C M be an open convex set and X a point-set vector fields on M . X is 
said to be monotone on Q, if 

(u^xpp 1 q) < (v^expy 1 p), p,q£&, u G X(p), v G X(q). 

Remark 2.1. This last definition has appeared in but also it is worth to point out that an 
equivalent definition has appeared in [35]. 

Let C C M be a set nonempty, convex and closed. It is well-known (see Corollary 3.1 of [1]) 
that for each p € M there exists a unique element p € C such that 

(expr 1 p, expr 1 q) < 0, q G C. 
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In this case, p is the projection of p onto the set C which we will denote by Tlc(p)- 

Let D C M n be a convex set, and p E D. Following [2], we define the normal cone to D at p by: 

iV D (p) := {w 6 T p M : (w, exp^ 1 q) < 0,q € D}. 

The previous definition holds just when M is of the Hadamard type. A more general definition has 
appeared in [25J . 

3 Proximal Point and Weak Sharp Minima on Riemannian Man- 
ifolds 

The definition of weak minima sharp as well as some of the basic related results that are used in this 
paper were introduced, in the Riemannian context, by Li et al. [25J for constrained optimization 
problems in the Riemannian context. Assuming that the solution set of Problem ([1]) is a set of weak 
sharp minimizers associated to the problem in question, we proved that the proximal point method 
([2]) has finite termination. We recall first some basic facts on the sequence generated by proximal 
point method ([2]) (see, for instance, [1]). In the remainder of this paper / : M — >• K represent a 
convex function, U denote the solution set of Problem [1] and {p k } is the sequence generated by ([2]). 
Moreover, we suppose that U is nonempty and closed in M. 

Assuming that the sequence {p k } is well defined, it follows that 

{\ k /2)d 2 {p k+ \p k )<f{p k )-f{p k+1 ). 

Hence, if / is bounded below, then 

J> fc /2)dV + V) < +°°- (4) 

fc=0 

The following proposition gathers some of the main results of pQ associated to the sequence {p k }. 
Proposition 3.1. If M is a Hadamard manifold, then the following statements hold: 
i) {p k } is well defined and is characterized by 

^k(expp^ +i p k ) € df(p k+1 ); (5) 

V Ysk=o = +°° an d th e solution set of the problem ([T]) is nonempty, then the sequence 
{p k } converges to a solution of the problem ([1]). 

Note that if A_, A + € M are such that < A_ < < A + , k G N, in particular Ylt=o ^Mfc = +°°- 
Assuming that < A_ < A^ < A + , k £ N and taking in consideration that || exp" 1 ^!! = d(q,p) (see 
Section f2.2j) . from combined with item i) of the previous proposition, it follows that 
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Corollary 3.1. If {(p k ,g k )} G TM (TM denote the tangent boundle, see Section^) is a sequence 
such that g k+1 = \k(exp~^ +1 p k ) G df(p k+ i), then 



+00 



^ll/H < +00. 



k=0 

In particular, for each a > there exists ko G N such that \\g k °\\ < a. 

Definition 3.1. The set U is said be the set of weak sharp minimizers for Problem{]\ with modulus 
a > if 

f(q)>f(p)+adu(q), p G U, q G M. 

Next proposition establish one characterization of the set of weak sharp minima on Riemannian 
manifolds. 

Proposition 3.2. A necessary and sufficient condition for U be the set of weak sharp minima for 
Problem [7] with modulus a > is that 

aB p nNu(p)cdf(p), P GU. 

Proof. Note that Nm(p) = {0} and df{p) is a closed set, for each p G M. Hence, the closure of the 
set df{p) + Nm{p) is equal to df(p), and the proof of this proposition follows from Theorem 4.6 
and 5.5 of EH. □ 



The following proposition is the key of our paper, since it is fundamental in the proof of the 
central result of this paper. 

Proposition 3.3. Suppose that U is the set of weak sharp minima for Problem U\ with modulus 
a > 0. If there exist q G M and w G T q M such that \\w\\ < a and w G df(q), then q G U. 

Proof. Let us suppose, by contraction, that q £ U and define 

exp^ 1 q 

Q = a p , p:=Uu(q). 
d{q,p) 

From the definition of p, it follows that 

(expp 1 q, exp~ 1 q) < 0, qeU. 

Combining definition of q with last inequality and definition of the normal cone Njj(p), we obtain 
q G Njj(p). Hence, since q G aBq, we conclude that 

q G aB p r\Njj{p). 
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Now, using that U is the set of weak sharp minima for Problem [T] with modulus a > 0, p G U and 
Proposition 13.21 last inclusion implies that q G df{p). Recall that df is a monotone field. So, since 
w G df(q) and g G df(p), we have 

(?, exp p 1 g) < (w, - expg 1 p). (6) 

Taking account that (w, — exp" 1 p) < \\w\\ || exp" 1 p\\, || exp" 1 ^ = d{q,p) and definition of q, from 
inequality ©, we obtain 

\\w\\ > a, 

which contradicts the inequality \\w\\ < a, and the desired result is proved. □ 

Next we present the central result of this paper. 

Theorem 3.1. Suppose that U is the set of weak sharp minima for Problem^ with modulus a > 
and let p° £ W 1 . If is a sequence of real numbers and A_, A + positive constants such that 
A- < Afc < A+, k G N, then the proximal point method terminates in a finite number of iterations. 

Proof. The proof follows immediately from Proposition l3 . 1 1 combining with Corollarv l3.1l and Propo- 
sition 13.31 □ 



4 Final Remarks 

In this paper we recall the notion of weak sharp minima for unconstrained optimization problem on 
Riemannian manifolds and we explored properties of weak sharp minimum on Hadamard manifold 
to establish the finite termination of the proximal point method. 
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